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Dynamics of hot gauge theories* 
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A brief overview is given of recent progress in understanding the dynamics of hot gauge theories. 
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1. Introduction 

Much progress has been made on understand- 
ing static equilibrium properties of gauge theories 
at finite temperature. This includes perturbative 
results valid at asymptotically high temperature 
as well as accurate results on phase struc- 
ture, thermodynamics, correlation lengths and 
other observables from numerical simulations 5 . 

Far less progress has been made on dynamic 
properties of hot gauge theories, such as equi- 
libration rates and transport properties, despite 
the fact that these properties are of direct inter- 
est in applications to both heavy ion collisions 
and early universe cosmology. The reason is sim- 
ple. Static properties may be extracted directly 
from the Euclidean theory, to which the whole 
panoply of modern theoretical tools (numerical 
simulations, renormalization group methods, ef- 
fective field theories, ...) may be applied. Dy- 
namic properties require analytic continuation of 
thermal correlation functions back to Minkowski 
space, or equivalently a functional integral formu- 
lation with a complex measure involving a non- 
trivial contour in the complex time plane |fi!7| . 

Nevertheless, considerable progress has been 
made in recent years in understanding dynamic 
processes in very hot, weakly coupled gauge theo- 
ries!]] The utility of theoretical results for asymp- 
totically high temperature to data obtained at 
current or future heavy ion experiments is an 
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quarks. For electroweak theory, this means T . 



open question. But at the very least, the asymp- 
totically high temperature regime is an instruc- 
tive theoretical laboratory which serves as a 
warm-up for efforts to understand more realistic 
non-equilibrium situations. 

In the asymptotically high temperature regime 
where the running coupling g(T) is small, one 
may disentangle a variety of phenomena which 
depend on different characteristic spatial or tem- 
poral scales. Dynamic phenomena display a 
richer set of characteristic scales than do static 
equilibrium properties. In a non-Abelian gauge 
theory like QCD, these include: 

T _1 The energy of a typical quark or gluon is 
0(T). Hence, T _1 sets the scale for the de 
Broglie wavelength of typical excitations. 

(gT)" 1 Electric fields are Debye screened on this 
length scale. The corresponding frequency, 
O(gT), is the scale of thermal corrections to 
the energy of typical excitations, as well as 
the plasma oscillation (plasmon) frequency. 

(g 2 T lnl/g) -1 This is the characteristic "color 
coherence length" — the maximum length 
over which a quark or gluon can be regarded 
as having a definite color. The factor of 1/g 
inside the logarithm comes from a ratio of 
the (gT)' 1 Debye length and the {g 2 ^ 1 
magnetic length. 

(g 2 T)~ 1 The amplitudes of magnetic field fluctu- 
ations on this length scale (or longer) are 
sufficiently large that their dynamics be- 
comes non-perturbative. 

(g 4 T In 1 1 g)^ 1 This is the characteristic large an- 
gle scattering time for quarks or gluons — 
the mean time for their direction of motion 
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to change by 0(1). This is also the charac- 
teristic relaxation time for non-perturbative 
magnetic fluctuations whose wavelengths 
are 0[(g 2 T)^}. 

For sufficiently small coupling g, the mean free 
path between scattering events which change a 
quark or gluon's color is longer than their de 
Broglie wavelength by a parametrically large fac- 
tor of 1 / g 2 In g~ x , and the mean free path between 
large angle scattering events is longer still by an- 
other factor of l/<? 2 . Consequently one should 
regard quarks and gluons with typical "hard" mo- 
menta (p ~ T) as well-defined, weakly interacting 
quasiparticles. 

Since the phase space of hard excitations is 
parametrically large compared to that of soft 
(p ~ gT) excitations, most physical observ- 
ables are predominately sensitive to the prop- 
erties of hard excitations. This includes bulk 
thermodynamic quantities like pressure or energy 
density, the Debye screening length, equilibra- 
tion rates and transport coefficients, the photo- 
emission rate, and many others. But there are 
important exceptions. In the high temperature 
phase of electroweak theory, for example, the rate 
of baryon number violation [Hj , and the wall ve- 
locity of a bubble nucleated at a first order phase 
transition depend on the low- frequency dy- 
namics of non-perturbative magnetic fields. 

2. Quasiparticle scattering 

Understanding the various scattering processes 
which can affect a hard quasiparticle is a prereq- 
uisite for calculating any observable which probes 
quasiparticle dynamics. These scattering pro- 
cesses include the basic 2 «-> 2 particle processes 
of gauge boson exchange, Compton scattering, 
and pair annihilation illustrated in Fig. ^ For a 
hard scattering with momentum transfer q ~ T, 
the explicit g 2 in the amplitude, plus dimensional 
analysis, implies that the rate is of order g A T. 
But a soft scattering, with momentum transfer 
q ~ gT, has a much larger 0(g 2 T) rate due to the 
1/q 2 behavior of the exchanged gluon propagator, 
which reflects the long-range nature of Coulomb 
interactions. This enhancement of soft scattering 
is essentially cut off by Debye screening at the 




Figure 1. 2 <-> 2 particle scattering processes. 




Figure 2. Near-collinear bremsstrahlung and in- 
clastic annihilation processes. 



gT scale. Note that a soft scattering event causes 
only a small 0(1/ g) change in a quasiparticle's 
direction, but will make an O(l) change in its 
color. 

The 0(g 4 T\ng~ 1 ) large angle scattering rate 
is the same as the two-body hard scattering rate 
except for the Inl/g factor. This arises because 
an 0(1) change in angle need not be produced by 
a single hard scattering — it may also occur via a 
sequence of many soft scatterings each making a 
small change in the particle's direction, which add 
incoherently to produce an 0(1) deflection. This 
gives rise to a log enhancement which is cut off at 
the soft gT scale. And the 0[{g 2 T In l/g) _1 ] color 
coherence length is just the inverse of the two- 
body soft scattering rate again up to a log factor 
that arises from logarithmic sensitivity to trans- 
verse magnetic scattering with q <C gT, which is 
only cut off by non-perturbative physics |1()| . 

In addition to these 2 <-» 2 particle processes, it 
is also important to consider the bremsstrahlung 
and inelastic annihilation processes illustrated in 
Fig- lTl413| . These amplitudes contain an ex- 
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plicit factor of g 3 , so on would naively expect an 
0(g e T) rate which would be suppressed relative 
to the g T hard scattering rate. However, if one 
examines the near-collinear, soft-exchange region 
in which the momentum transfer q ~ gT and the 
angle cj> ~ g, then one finds that these rates have 
a l/g 4 T 4 enhancement from the exchanged soft 
gluon propagator, a further 1/g T 2 enhancement 
from the nearly on-shell internal quark propaga- 
tor, and a g 2 suppression from the near-collinear 
kinematics at the photon vertex. When combined 
with the explicit g e factors plus a g 4 T 4 phase 
space suppression from the restricted kinematics, 
one finds an 0(g 4 T) rate — exactly the same as 
the hard scattering rate. 

Consequently, a hard quark (or gluon) moving 
through the plasma can "fission" into a nearly 
collinear pair of hard excitations, or "fuse" with 
a another nearly collinear hard excitation, at a 
rate which is parametrically the same as the two 
body scattering rate. These processes may be 
regarded as 1 «-> 2 particle processes which are 
normally forbidden by kinematics, but which be- 
come kinematically allowed when accompanied by 
a soft transfer of momentum to other particles 
in the system. However, if one soft scattering 
with the rest of the system can take place dur- 
ing a near-collinear bremsstrahlung or annihila- 
tion process, so can two or more. The l/g 4 T 2 
enhancement from the internal quark propaga- 
tor is a sign that virtual intermediate states in 
these processes are living for a parametrically 
long time of order \jg 2 T — which is not small 
compared to the color coherence time. Therefore, 
multiple soft scatterings occurring during these 
near-collinear processes cannot be neglected and 
will produce an O(l) suppression in the resulting 
rate. This is known as the Landau-Pomeranchuk- 
Migdal (LPM) effect QUE]. 

The net result is that a leading- order calcula- 
tion of any observable which is sensitive to the 
dynamics of hard quasiparticles must correctly 
incorporate both the appropriate two-body scat- 
tering processes as well as LPM-suppressed near- 
collinear emission processes. 
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Figure 3. Total constant under the log, C(k/T), 
for photon emission from hot QCD. (Reproduced 
from Ref. |2Z|.) 



3. Photon emission rate 

The photon emission rate from hot QCD is one 
quantity which has received considerable theoret- 
ical attention |llH13ll8l42fi| . due in part to the 
hope that it will provide a useful probe in heavy 
ion collisions. A complete leading-order evalu- 
ation of this emission rate (in equilibrium) has 
only recently been completed [23 . The differen- 
tial emission rate for photons of momentum k has 
the form 

^ = B(k) pnCT/moc) + C(k/T) + 0(g s )} , (1) 

where B{k) = ~ oem ct s T 2 (J2 s q 2 ) k/{e k / T +l) , 
q s is the charge assignment of quark species s, 
and m 2 ^ = g 2 T 2 /3 is the asymptotic thermal 
quark mass. The function C(k/T), shown in 
Fig- El is the "constant under the log" ; it is a non- 
trivial function of k /T but it is independent of the 
strong coupling g s , whereas hi(T /moo) ~ m (l/<7s) 
since = 0(g s T). The lowest order two body 
processes qq — » gj and qg — > (77 generate the 
hxiT/nioo) term and part of C(k/T); these con- 
tributions were evaluated in |18I19| . The near- 
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collinear processes only affect the non-logarithmic 
term C(k/T). Their contribution to the leading- 
order emission rate was recently shown to in- 
volve the solution of a non-trivial integral equa- 
tion which incorporates the effects of multiple soft 
scatterings during the emission process, dynami- 
cal screening in the plasma, and thermal correc- 
tions to the quasiparticle dispersion relations [2S| • 

4. Transport coefficients 

Transport coefficients, such as shear viscosity 
or flavor diffusion constants, are physical observ- 
ables of obvious interest which are sensitive to the 
dynamics of quasiparticles. Specifically, they de- 
pend on the rates at which non-equilibrium fluc- 
tuations in the phase space distribution of quasi- 
particles relax back toward equilibrium. And 
these rates depend on exactly the same scattering 
processes depicted in Figs. ^ and [21 Parametri- 
cally, these transport coefficients are proportional 
to the large angle scattering time and hence scale 
as (g 4 In l/g) -1 . The fact that transport coef- 
ficients grow linearly with mean free time is well 
known from the simple Drude model of conductiv- 
ity; the current induced by an applied field is pro- 
portional to the length of time between collisions 
during which charges are accelerated by the field, 
ft is the large angle scattering rate which is rele- 
vant, not the much faster rate of soft scattering, 
because these transport coefficients all involve the 
flux of gauge invariant conserved densities (i.e., 
momentum density, isospin density, etc.), and the 
contribution of a quark or gluon to these fluxes 
is not significantly changed if it undergoes a scat- 
tering with tiny momentum transfer □ 

These transport coefficients may be formally 
defined by Kubo formulas relating them to the 
zero frequency limit of a current-current or stress 
tensor-stress tensor Wightman correlation func- 
tion. The inverse dependence of the coupling 
constant is an obvious sign than an infinite num- 
ber of Feynman diagrams must contribute to the 
leading order result |f| However, the complicated 



2 This is not true for the "color conductivity" of a non- 
Abelian plasma, which is sensitive to very soft scattering 
processes. 

3 The usual loop expansion breaks down because one is 



set of diagrams which contribute to the leading- 
order result may be summed up by a suitable 
integral equation, which is precisely a linearized 
Boltzmann equation (projected onto a particular 
symmetry channel) for small perturbations in the 
quasiparticle phase space densities. 

As a specific example, the flavor diffusion con- 
stant characterizing the relaxation of fluctuations 
in isospin or strangeness density (in the high tem- 
perature limit where quark masses are negligible) 
has the form 



D = A/{g*T [Ins- 1 + 0(1)]} 



(2) 



where A is a pure number (depending on the 
number of active quark flavors). Just as for the 
photon emission rate, the leading-log coefficient 
A is determined by 2 <-> 2 scattering processes, 
while the 0(1) constant "under" the log also de- 
pends on near-collinear emission and absorption 
processes in the presence of multiple soft scat- 
terings. Quite a few efforts have been made to 
evaluate QCD transport coefficients like D in a 
"leading-log" approximation (where one pretends 
that lnl/(7 s 3> 1) by solving a linearized Boltz- 
mann equation incorporating appropriate 2 <-> 2 
particle scattering rates |29H35| . Nevertheless, 
until recently almost all reported results were in- 
correct due to a failure to appreciate that Comp- 
ton scattering and pair annihilation processes, 
in addition to gluon exchange, contribute to the 
leading- log result. This has been remedied in re- 
cent work |8fi| which found 
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(3) 



24 + 4N f + 7T 2 ' 

for QCD with Ni flavors^ The three terms in 
the denominator arise, in order, from t-channel 
gluon exchange between a quark and a gluon, t- 
channcl gluon exchange between two quarks, and 
Compton scattering or annihilation to gluons. 

Analogous leading-log results for shear viscos- 
ity and electrical conductivity may also be found 
in [H]. 

evaluating the correlation function at an exceptional point 
in momentum space, namely uj, k — * 0. 
4 This is an approximate form, accurate to within a frac- 
tion of a percent, which results from using a one-term 
variational ansatz in the relevant integral equation 36 . 
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Of course, a leading-log result, in which the 
constant under the log is totally undetermined 
(and hence relative corrections suppressed only 
by l/lng" 1 are neglected) is unlikely to provide 
a useful prediction in any realistic theory. At the 
very least, one would like to obtain a complete 
leading-order result, in which neglected effects are 
suppressed by at least a factor of the coupling g s . 
Doing so requires the correct inclusion of near- 
collinear gluon emission or absorption processes 
analogous to the near-collinear processes relevant 
for photon emission. It should be possible to 
augment the linearized Boltzmann equation with 
effective 1 <-> 2 particle scattering terms which 
correctly reproduce these near-collinear reactions, 
and whose transition rates incorporate the correct 
LPM suppression effects. (However, this will re- 
quire solving a non-trivial integral equation just 
to determine the kernel to be used in another inte- 
gral equation.) Such work is currently in progress. 

5. Soft gauge field dynamics 

As noted earlier, a few important observables 
like the rate of baryon number violating transi- 
tions (also known as the "topological transition 
rate") are not primarily sensitive to the dynam- 
ics of hard quasiparticles, but instead probe low 
frequency, long wavelength gauge field dynamics. 
The starting point for understanding this regime 
is the observation that the relevant degrees of 
freedom, namely modes of the gauge field with 
k, u> -C T, will have parametrically large occupa- 
tion numbers since the Bose distribution UbipS) ~ 
T/ui for u> CT. Consequently, these modes may 
be viewed as classical fields [37|. But these soft 
modes of the gauge field are driven by the color 
current generated by all the hard quasiparticles 
in the plasma. If one splits the theory into hard 
and soft degrees of freedom, one may formulate a 
Boltzmann- Vlasov kinetic theory which describes 
the propagation of ultrarelativistic hard excita- 
tions in the background of a long wavelength clas- 
sical gauge field, together with the non-Abelian 
Maxwell equation D v F tJ,l/ = characterizing the 
reaction of the hard quasiparticles back on the 
soft gauge field. This kinetic theory (linearized in 
the deviation of quasiparticle distributions away 



from equilibrium) is a formulation of the well- 
known hard thermal loop (HTL) effective the- 
ory I2BHHI- 

As shown by Bodeker |41j , one may systemati- 
cally integrate out the effects of fluctuations down 
to a scale fi <C gT, and derive an effective theory 
for the soft gauge field alone. Of course, since 
hard quasiparticles propagate as nearly free par- 
ticles, with definite color, over distances up to the 
color coherence length 7 -1 = 0[(g 2 Tlnl/g) ], 
this effective theory will be non-local on the scale 
of 7 _1 . But if one restricts attention to distances 
large compared to the color coherence length, 
or wavenumbers k <C 7, then one may formu- 
late a consistent local effective theory. The re- 
sult is a stochastic theory which simply combines 
Ampere's law, Ohm's law, and the fluctuation- 
dissipation theorem [41142) . 



DxB = ffE-(. 



(4) 



The single parameter a is the color conductiv- 
ity, and C is Gaussian noise with a variance 

As noted in the introduction, the color coher- 
ence length 7 _1 is parametrically smaller than 
the non-perturbative magnetic length by one fac- 
tor of 1/lng . Consequently, if one applies 
this effective theory to non-perturbative g 2 T scale 
physics, then corrections to this effective theory 
will be suppressed by powers of 1/lng . In 
fact, it is possible to show 03] that corrections 
are suppressed by two powers of 1/lng -1 pro- 
vided one determines the correct value of the color 
conductivity by appropriately matching the long 
distance effective theory to the underlying HTL 
theory which in turn is derived straight from hot 
QCD. One finds (for an S\J(N C ) theory) 



3iV c a T 



In 



rnp 
70) 



C 



(5) 



where mo is the leading-order Debye mass (in- 
verse screening length), = N c aThi(rnu/ A*)j 
[i is a renormalization scale which should be cho- 
sen to be of order 7, and C = 3.041. The leading- 
log coefficient was found in [41144) , while the con- 
stant under the log, C, is the result of the very 
tricky matching calculation in [45). 
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In the effective theory (@J pure dimensional 
analysis shows that the the topological transition 
rate (per unit volume) has the form 

r = k {aTf/a (6) 

with k a dimensionless pure number which de- 
pends on the non-perturbative dynamics of the 
theory. But this effective theory is a super- 
renormalizable, UV finite theory which may be 
discrctized on a spatial lattice and numerically 
simulated in a completely unambiguous fashion. 
(In fact, in A n = gauge, the effective theory 
(@J is precisely the stochastic quantization of 3d 
Yang Mills theory.) Such a numerical simula- 
tion was performed in |46| . with a result that 
k = 22.6 ±1.5. 

The topological transition rate has also been 
extracted from real time simulations 47 48 of 
two different more microscopic (i.e., less "effec- 
tive") theories which are lattice versions of the 
Boltzmann-Vlasov theory described earlier. Sur- 
prising good agreement was found between the 
results obtained using the long distance effective 
theory |0J and both of these more microscopic 
formulations. As a result, the non-perturbative 
high temperature baryon violation rate is now a 
satisfyingly well determined quantity |48j . 

6. Open questions 

Many aspects of hot gauge field dynamics offer 
opportunities for further progress. Examples of 
problems which can be addressed using perturba- 
tive methods include the following. 

1. Perform complete leading-order evaluations 
of QCD transport coefficients such as shear 
viscosity, flavor diffusion, and electrical con- 
ductivity. 

2. Evaluate bulk viscosity in hot QCD. Bulk 
viscosity does not receive a leading-log con- 
tribution from 2^2 particle processes, and 
at present there are no published results de- 
riving even the parametric dependence on 
coupling in a gauge theory. 

3. Compute any transport coefficient beyond 
leading order, even in a non-gauge theory. 



Can a valid effective theory still take the 
form of kinetic theory? 

Other natural topics involve the development 
or improvement of methods for studying real time 
dynamics outside of perturbation theory. To list 
just a few goals: 

1. Test the domain of utility of leading order, 
or leading-log perturbative results for ob- 
servables like transport coefficients. This 
has been accomplished, so far, only in one 
special case, namely the iVf — ► oo limit |49| . 
where the next-to-leading log approxima- 
tion works surprisingly well as long as the 
coupling g 2 N[ is not so strong that its scale 
dependence becomes a dominant effect. 

2. Develop a practical scheme for extracting 
transport coefficients from Euclidean lattice 
gauge theory simulations. This would in- 
volve extracting an estimate for the zero 
frequency slope of a spectral density from 
knowledge of a Euclidean correlator. 

3. Extract more physics from real-time clas- 
sical lattice gauge theory simulations. It 
should be possible to exhibit the presence of 
over-damped low frequency gauge field dy- 
namics in observables other than the topo- 
logical transition rate. Current efforts in 
this direction have found somewhat per- 
plexing results [5T3| . 

4. And last but not least, develop better meth- 
ods for exploring the real time dynamics 
of systems with substantial departures from 
equilibrium. A couple of recent steps in this 
direction include |51l52j . 
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